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Abstract. - Charge, spin, and orbital degrees of freedom underlie the physics of transition metal 
compounds. Much work has revealed quantum critical points associated with spin and charge 
degrees of freedom in many of these systems. Here we illustrate that the simplest models that 
embody the orbital degrees of freedom - the two- and three-dimensional quantum orbital compass 
models - exhibit an exact quantum critical behavior on diluted square and cubic lattices (with 
doping 5 = 1/4 and 5 = 1/2 respectively). This raises the possibility of quantum critical points 
triggered by the degradation of orbital order upon doping (or applying pressure to) such transition 
metal systems. We prove the existence of an orbital spin glass in several related systems in which 
the orbital couplings are made non-uniform. Moreover, a new orbital Larmor precession (i.e., a 
periodic change in the orbital state) is predicted when uniaxial pressure is applied. 



Introduction. The interplay between charge, spin, and 
orbital degrees of freedom is a key ingredient underlying 
the physics of transition metal compounds. The elec- 
tronic orbital degrees of freedom often allow for coop- 
erative effects leading to well defined spatial orderings. 
Generally, crystal field effects split the five n = 3 d- 
wave orbital wavefunctions into a triplet (i 2g orbitals of 
the xy ee (y 2i _ 2 - F 2>2 )/V2, yz = (Y 2 ,-i + Y 2 ,i)/V2 
or xz ee (i2,—i — l2,i)/v2 types) and an e g doublet 
(spanned by orbitals whose angular dependence is of the 
oz 2 - r 2 ee Y 2 , Q and x 2 - y 2 = (F 2 ,_ 2 + Y a>2 )/V2 forms). 
Colossal magnetoresistance (CMR) - the sharp decrease 
of resistivity with magnetic field - occurs in some of the 
compounds that display orbital orders, e.g. the manganite 
LaMn03 [1] and its derivatives. In LaMn03, the orbital 
order is of a simple character: the state of the single out- 
ermost electron in one ion may be that of, say, the n = 3 
d-wave state of the 3z 2 — r 2 type while it is 3a; 2 — r 2 on 
a neighboring ion and so on in a staggered fashion within 
each plane. Other prominent systems that exhibit various 
types of orbital order are found among the vanadates (e.g., 
V 2 3 [2], LiV0 2 [3] and LaV0 3 ) [4] and cuprates (e.g. 
KCUF3) [5]. Orbital ordering can be observed via orbital- 
related magnetism and lattice distortions or by resonant 
X-ray scattering techniques [6] . 

Although there are numerous works on quantum crit- 
icality in systems with various electronic phases [7, 8] 



there is little prior work [9-12] on systems in which the 
orbital ordering temperature veered to zero. We wish 
to motivate the following question: If quantum criti- 
cal points are associated with the degradation of mag- 
netic/charge/superconducting orders in numerous systems 
why can these not appear in orbitally ordered systems? 
Such quantum critical points are associated with the 
degradation of orbital order. This path has not been 
followed despite the success in finding quantum critical 
points and novel low-T transitions in spin and other elec- 
tronic systems. We demonstrate the existence of exact 
quantum criticality driven by orbital order in the sim- 
plest of all orbital only models: the D = 2 and 3 dimen- 
sional quantum Orbital Compass Models (OCMs). These 
systems rigorously exhibit order in their classical limit 
[13-15]. This order is expected to be fortified by quantum 
effects. Indeed, numerically orbital order was detected in 
the undiluted D = 2 OCM [11]. We will show that, at a 
prescribed doping (dilution) of magnitude 5 = 1/4, this 
system displays quantum critical correlations. This quan- 
tum critical phase may be driven away by applying pres- 
sure/strain (i.e. varying the size of the exchange constants 
corresponding to different spatial directions). Similarly, 
we show that a particular set of anisotropic couplings of 
the D = 3 OCM leads to quantum criticality on a lattice of 
doping S = 1/2. In the particular case that the orbital ex- 
change couplings become random (i.e., non-uniform), we 
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demonstrate that an orbital spin glass may be generated. 
Furthermore, for any orbital system, when uniaxial pres- 
sure is applied the phenomenon of orbital Larmor preces- 
sion develops, i.e., the orbital states change periodically 
in time as a consequence of the applied pressure. 

Orbital only models. The orbitals in the transition 
metal oxides are modeled by S = 1/2 SU{2) degrees of 
freedom [1,16]. These pseudo-spins appear in the orbital 
only Jahn- Teller (JT) interactions and spin-orbital Kugcl- 
Khomskii (KK) Hamiltonians depicting exchange. The 
simplest Hamiltonian linking interactions with the direc- 
tionality of the orbitals in space is the OCM [16]. The 
S — 1/2 anisotropic D = 2,3 anisotropic OCMs [17] on 
^Vocm = Yifi sites of a square or cubic lattices respec- 
tively are given, in term of Pauli matrices with \x = x,y 
(D = 2) or fi = x,y,z (for D = 3), by 

ffocM = -£.vK+, ti . (1) 

3 

Interactions involving the cc-component of the pseudo-spin 
Sj = o~j /2 occur only along the spatial x direction of 
the lattice. Similarly spatial direction-dependent spin- 
exchange interactions appear for the y (or z) components 
of the spin. In D = 3, this model emulates the direc- 
tional character of the orbital related interactions trig- 
gered by JT effects and the orbital only component of the 
KK Hamiltonian [1, 13, 14, 16] in materials with a single 
electron in a t^ g level. Varying the strengths of the cou- 
plings {J^,} relative to one another emulates the effects 
of uniaxial strain and pressures on the orbital-only com- 
ponent of the orbital-dependent spin-exchange (the KK 
Hamiltonian). Similarly, the influence of uniaxial pres- 
sure on the JT interactions may be mimicked by adding a 
pseudo-spin magnetic field linearly coupled to the orbital 
pseudo-spins. For instance, uniaxial pressure along the 
cubic z axis favors orbital states that have little extent 
along the z axis. In the case of e g orbitals, such a uniaxial 
pressure favors the planar d x 2_ y 2 states (represented by 
a z = —1 [1, 13, 14, 16]) vis a vis the cigar-shaped d 3z 2_ r 2 
state (depicted by a z = 1), that is elongated along the z 
axis. 

Detailed calculations were performed on orbital systems 
such as the OCM [18]. We will mainly focus on the 
D = 2 variant of the OCM. By the mappings of [17], the 
D = 2 OCM is dual to (p + ip) superconducting Joseph- 
son junction arrays on a square lattice. Recent theoretical 
work suggests similar interactions in p-band Mott insula- 
tors [21]. 

The planar orbital compass model. The D = 2 
OCM has the following symmetries [19, 20] 6^ = 
Iljec ' f° r l 1 = x iV with C M _L e M axis. On a 

torus, these operators arc defined along toric cycles. As 
O m involves ©(L 1 ) sites, they constitute d = 1 symme- 
tries [19,20]. Such symmetries cannot be broken at finite 
temperatures [20]. These are also symmetries of the di- 
luted model that we will describe below. At the isotropic 



point J x = J y it has an additional global reflection sym- 
metry 

Reflection = JJ e*^ W+"i"> , (2) 
3 

which may be broken at finite temperatures [22]. 

As we will prove below, the diluted OCM exhibits a 
quantum critical point at the very same location ( J x = J y ) 
where a finite temperature transition is expected to oc- 
cur in the undiluted OCM. It is natural to anticipate 
that the ordering temperature drops monotonically with 
doping [11] with, as wc will show, the system becoming 
quantum critical for a prescribed dilution of S = 1/4 
(see Fig. 1). Away from the isotropic point J x = J y , 
all symmetries are of the O m type. At finite tempera- 
tures, these symmetries cannot be broken and only topo- 
logical quantum order may be possible [19,20]. The al- 
gebra satisfied by the bond operators A^ = cfoj 1 , for 
pairs of sites in which j — i = e M , is particularly simple: 
{A^AD = {A^Al,} = 0, (j ± k + i), [A^A^] = 0, 

(A^ = 1, [A^A^] = 0, (/x £ J , i,j ± k,l). Apart 
from these algebraic relations, there are no additional con- 
straints that link the bonds {^j}- 

A doped orbital compass model. Let us now consider 
the Hamiltonian of Eq. (1) on the diluted lattice of Fig. 
1 (lower panel) with N — I^ocm sites, and call this new 
Hamiltonian -Hdocm- The lattice of Fig. 1 corresponds 
to a doped (or diluted) system in which 5 = 1/4 of the 
sites have been replaced by an inert site (dopant). After 
doping, in addition to the symmetries of O m and Reflection, 
a new gauge symmetry emerges 

6 a = C7 X a , [H DO cM,d a } = 0, (3) 

where a represents sites which are two-fold coordinated, 
i.e. z = 2 (denoted by squares in Fig. 1). Hence- 
forth we will denote the set of z = 2 sites by r^. This 
gauge symmetry allows to decompose the total Hilbert 
space into orthogonal subspaces of equal dimensionality 

^ = (Bt=i Me, where dimJQ = 2 2Ar / 3 . In each sec- 
tor for all states \4> n )e S ^e, we have that a^\4> n )i = 
Va\<f>n)i , Va = ±1- There are 

2 N/3 

sequences of eigen- 
values of these operators. Thus, each subspacc is labeled 
by a particular string of ±1, i.e. !Ke = J{, ,»/3, with for 
example {n a } = +, +, +,—,+,—,•■■ , +, +. Each of the 
Hilbert subspaces %i spans the (2N/3) spins on the re- 
maining 2iV/3 (z — 3) sites. These sites lie along vertical 
columns separated by one another by intervening columns 
of z = 2 sites (see Fig. 1). Let us label the set of all 
three-fold coordinates sites b by f2 3 (the sets and 
are disjoint fi 2 H ^3 = {0})- The algebra of the bonds 
{A^j} in each of the projected subspaces 'Kt is unchanged 
relative to that defined on the full Hilbert space. To prove 
this, let us define the projection operators 

JV/3 

Pi = tl i^P") ' P * = Pe > ^CM,Pi] = 0. (4) 

a=l 
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Fig. 1: Doped D — 2 orbital compass model. On each vertex 
of the graph there is an S — 1/2 degree of freedom indicated 
by a spin-1/2 operator Sj. The top panel represents the un- 
doped orbital compass model (z = 4) with circles indicating 
the dopant sites. Vertical links have strength J y while horizon- 
tal links have strength J x . The lower panel displays the result- 
ing doped model where the graph still remains two-dimensional 
(brick-wall topology) but with two types of vertices: one shares 
three neighbors (z = 3) while the other (squares) is only con- 
nected to two neighbors (z — 2). 



Within each of the projected subspaces, the bond oper- 
ators A^j = PiAjjPi satisfy relations identical to those 
above, i.e. {A^A v ik } = {A^A^} = 0, (j ? k jt i), 
[At v A^] = 0, P t (A^P t = u [A^Ai] =0,(n? 

fjf, i,j 7^ k,l). These relations follow directly from the 
fact that [Pi,A x tj \ = 0, since [P e , a]} = 0, and [P e , AVj] = 0, 
because all y-type bonds A\- have their support entirely 
in the domain of z = 3 sites i,j C fl^. Consequently, the 
bond algebra does not change. 

We now exactly solve the doped OCM (DOCM). The 
algebraic relations for the interaction terms A^ defined 
above are identical to those of the interaction terms in the 
D = 1 transverse field Ising model (TFIM). In both in- 
stances, the dimension of the representation of the algebra 



is the same as are all additional constraints between the 
interaction terms. In any given sector Jt^ (which amounts 
to a particular gauge fix), all bonds in the DOCM system 
are those of decoupled TFIM chains. The Hamiltonian of 
a transverse field (hj) Ising chain is 

"mv Xt''" 7 ; ' )• (5) 

3 

On the other hand, the DOCM Hamiltonian of Fig. 1 in 
the projected subspace, /^docm = PeHoocMPe, is 

i^DOCM =-^2 (jxVa of + J y cr^al +iy ) , (6) 
b 

where a denotes the z = 2 site a £ S^2 that is a nearest- 
neighbor of a given z — 3 site b G (a lies to the imme- 
diate left or right of b in Fig. 1). #docm represents the 
Hamiltonian of a set of D = 1 TFIM chains of length L 
correlated by the gauge fields r] a , but otherwise uncoupled. 
This is equivalent to say that -Hdocm represents a single 
D = 1 TFIM chain of length 2N/3 and transverse fields 
hb = JxVa- The spectrum of a TFIM chain is indepen- 
dent of the sign of hj . Given the unitary (and hermitian) 
operator Ur = Yijer > wnere T denotes an arbitrary set 
of lattice sites, UrHjf\M[{hj}]Ur = Hjf\M[{hj}], where 
hj = —hj if j G r, and hj = hj, otherwise. Thus, all 
different projected -Hdocm have exactly the same spec- 
trum (although their eigenstates are interchanged). As 
each sector 3-Q leads to an identical partition function, 
the resulting one for the DOCM is 

Z = tr M e-' 3ifD0 ™ =J2 trx t e-P Ro °™ = 2 N ^Z TFlM , (7) 

i 

where Ztfim is the partition function of the TFIM of Eq. 
(5) with 27V/3 sites. If the system is on a torus of size 
L x x L y we can impose periodic boundary conditions if 
the dimensions of the lattice are such that there is an 
even number of sites along a vertical column and a num- 
ber of sites that is a multiple of four along the horizontal 
direction, i.e. L x = 0(mod 4), L y = 0(mod 2). In each 
sector %i in addition to the algebraic relations defined 
above, the vertical bonds will satisfy the additional con- 
straint: Ub=i ^-b,b+e y = 1. Analogously, the TFIM of Eq. 
(5) in which the bonds satisfy the analogous constraint 
rib=i a b a b+e = 1 corresponds to a system of decoupled 
chains (formed by the z — 3 sites) in each of which there is 
a periodic boundary condition (the chain is of length L y ). 
For a system of decoupled L x /2 TFIM chains the par- 
tition function is Zjfim = {zjf\u) Lx ^ 2 , where ztfim(£k) 
is the partition function of a TFIM chain of length L y . 
i?TFiM can be diagonalized by a Jordan- Wigner transfor- 
mation followed by a Bogoliubov transformation, -Htfim = 
Sfe £ fc(7fe7fc — 1/2) , with 7 fc (7!) the fermionic annihila- 
tion(creation) operator for a fermion (Bogoliubov quasi- 
particle) of wavenumber k (the set of allowed wavenumbers 
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is {— 7r, — 7r + Ti/Ly, ■ ■ ■ , 0, • • • , it — ir/Ly}. The quasipar- 
ticle energies of a TFIM chain of length L y are given by 



€k = ZyJ'x + Jy ~ 2J x J y cosk, and its ground state en- 
ergy Eq = — £fc/2. Thus, the grand potential per site 



fi 2£ fcsTln2 2k B T ^ 

TV ~3i^ 3 317^ 1 



(8) 



with z = exp[/3£] the fugacity and £ the chemical potential. 
Thus, the DOCM of Fig. 1 has a quantum critical point 
at J x = ±J y (D = 2 Ising universality class). Using our 
mapping, in the DOCM, the correlator 



G(r,t) = (tT v (f,t)o*(0,0))x = (a y (f,t)a v (0, 0)) 



TFIM- 



(9) 



Due to the decoupling to TFIM chains, the last average 
is zero unless r lies along the y axis. With c = 2 J y (the 
lattice constant is set to one), for T — 0, we have at the 
quantum critical point G(r, 0) ~ (|y|/c) — 1 ' 4 5 a;) o, and is 
manifestly critical. At the quantum critical point for gen- 
eral T > temperatures, the Fourier transformed correla- 
tor is given by [7] 



G{k, to) cx 



Hi 
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(10) 



By our mapping, expressions are similarly available in the 
paramagnetic and ordered phases. The regular pattern of 
doped sites leads to the emergence of the gauge symmetry 
O q , crucial for the exact solvability of the model. However, 
essential to observe lower dimensional (d = 1) physics is 
the presence of the d = 1 O' 1 symmetries [19]. A modest 
amount of randomness in the doping (i.e. not with a regu- 
lar pattern) docs not invalidate the dimensional reduction 
argument, although we can no longer solve the problem 
exactly. This conclusion is supported by the numerical 
simulations of Rcf. [11]. 

Recently, non-uniform structures that may be generated 
by doping orbital systems were investigated in Ref. [23]. 

The doped D = 3 orbital compass model. Con- 
sider a cubic lattice in which consecutively layered planes 
(stacked along the cubic z direction) have the following 
form: ABCBABCBA ■ ■ ■ . Here, B planes are as in Fig. 
1. Planes of the A and C types are those in which the 
corresponding doped columns of plane B are fully doped. 
Furthermore, along the columns that are undoped in plane 
B, planes A and C are half-doped with every other site 
removed such that there is a relative shift, between the 
A and C planes, of one lattice constant along the y di- 
rection. The average dopant density on this lattice is 
5 = 1/2. Consider the D = 3 OCM of Eq. (1) on this 
lattice. The symmetry of Eq. (3) remains intact on all 
planes of the B type, while {<x^} for all undoped sites 
u in the A and C planes are local symmetries. Repli- 
cating the solution of the D = 2 case, we find that the 
D = 3 DOCM reduces to decoupled TFIM chains in the 



B planes as before with a transverse field of strength 
h = J J% + J| . This is so as along the undoped pla- 
nar columns of Fig. 1, each site now feels two transverse 
fields along the x and z directions. The additional inter- 
actions along the z direction originate from a site in the 
A or C planes. No additional x or y interactions appear 
within the A or C planes. The dispersion is now given 



by e k = 2^J Jl + J2 + J? - 2 J y yfjf+Ji cos k. This par- 
ticular lattice exhibits quantum critical points along the 
locus J x + = Jy. As a function of uniaxial pressure 
along the y axis, the point J y = J\[2 (with J x = J z = J) 
is quantum critical. For increasing/decreasing uniaxial 
pressures, and their influence on the orbital-dependent 
spin-exchange, G(r, t) exhibits the correlations of the or- 
dered/paramagnetic phases of the TFIM [24]. 

Uniaxial pressure on John-Teller orbital-only interac- 
tions: Orbital Larmor precession. Thus far, we exam- 
ined the (diluted) anisotropic OCM. The latter captures 
the effects of uniaxial strain and pressure as these apply 
to the orbital component of the orbital-dependent spin- 
exchange interactions. We now turn to the influence of 
the uniaxial strain/pressure on the direct (JT borne) or- 
bital interactions. Following our earlier discussion, for a 
uniaxial pressure/strain along the v space direction, these 
now generally lead to a Hamiltonian of the form 



OCM: strain 



OCM 



-Hp 



where the uniaxial pressure is represented by 



(11) 



(12) 



with no summation over v. For the spatial directions 
v = x or y (here one needs to perform a unitary trans- 
formation), the Hamiltonian of Eq. (11) can, once again, 
be exactly solved in the two- and three-dimensional dilu- 
tions considered hitherto. Let us consider the case v = x. 
In this case, the symmetries of Eq. (3) remain manifestly 
unchanged. Following our earlier derivation, we arrive at 
a simple generalization of Eq. (6) 



DOCM; strain 



Hi 



DOCM 



b 



Y.^^, (13) 



i.e., a transverse field Ising model with an effective trans- 
verse field hb = J x ila + P x while, as before, the Ising ex- 
change constant J = J y . Unlike Eq. (5), (13) does not 
exhibit a symmetry of the spectrum under local inversions 
of the type rj a — > — rj a . The general case of uniaxial strain 
cannot be exactly solved and leads to a richer system. 
We have seen that the pressure induced orbital-dependent 
spin-exchange interactions may, on their own, lead to a 
uniaxial pressure tuned quantum critical point. We can- 
not trivially solve for the thermodynamics once the direct 
orbital JT interactions of the form of Eqs. (12) and (6) are 
included. The crux is the trace over all {n a } that should 
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be performed which in effect sums over the partition func- 
tion of all subsystems with different transverse fields hi,. 
Hamiltonian (13) can be exactly solved within each {n a } 
sector. It may well be that the quantum critical behav- 
ior persists or that, e.g., a Griffiths-type phase appears in 
place [7]. 

We now comment on a related feature of the diluted 
OCM when spatially non-uniform couplings are present. 
Consider a general system in which the (random) exchange 
constants {J^j} vary from bond to bond (defined by (j, j+ 
in)). Replicating our derivation above, we find that after 
applying a uniform external field that emulates the effect 
of uniaxial strain/pressure, we arrive at a transverse field 
random Ising model (TFRIM) within each sector of {rj a }, 

^TFRIM = - (i J *.bVa + P x ] + Jy.b ^X+eJ ■ ( 14 ) 
b 

Iii Eq. (14), P x plays the role of a uniform transverse field 
while {J y ,j} are random longitudinal exchange constants. 
As the TFRIM exhibits spin glass behavior with well 
known characteristics, so may does our system. We may 
thus predict the appearance of an orbital spin glass [26]. 

A fundamental consequence of the form of the uniax- 
ial pressure/strain induced JT interaction is the appear- 
ance of orbital Larmor precession. This precession hap- 
pens even in the simplest bare orbital system. If only the 
consequences of uniaxial pressures along a direction v are 
considered, then for any orbital system, the Hamiltonian 
Hp of Eq. (12) leads to orbital Larmor precessions. That 
is, the orbital states of the electrons change periodically 
in time due to the application of pressure 

^- =1 axP, (15) 

with P = ^2 U P V e„, and 7 the effective gyromagnetic 
ratio of the orbital system. The precession times scale as 
|P| _1 . In principle, our new orbital precessions should be 
observable. 

Possible realizations of orbital ordered triggered quan- 
tum criticality. Our work shows that in orbital sys- 
tems, upon the application of external pressure (uniaxial 
strain) and varying doping, we may attain quantum crit- 
icality associated with the degradation of orbital order. 
What physical systems may display new orbital fluctua- 
tion driven quantum critical points? The Luttinger liquid, 
a D = 1 electronic liquid, and some variants thereof are 
prototypical examples of critical systems. (Indeed, our 
exact solution relied on a dimensional reduction to decou- 
pled chains ensured by the symmetries of these models.) 
In some physically probed orbital systems, e.g. [25], an 
effective reduction in dimensionality to precisely such a 
system occurs due to the arrangement of orbitals. KCUF3 
has a weakly coupled spin chain structure and offers an 
example of a system with aD = l (Luttinger liquid) like 
behavior observed over a broad range of temperatures, mo- 
menta, and frequencies [27]. The orbital only component 



of the orbital dependent spin-exchange shares the same 
features and anisotropic form as the orbital only interac- 
tions [1,5]. The directional character of the orbitals in 
KC11F3 triggers a quasi one-dimensional behavior of the 
spins. Ionic substitution of KCUF3 (to KCui-sTingF^) [9] 
revealed a vanishing of the orbital ordering temperature 
for a doping of S = 1/2. A reduction in the orbital order- 
ing temperature occurs as Zn does not have an orbital de- 
gree of freedom. References [10,11] noted that this doping 
fraction is smaller than the one needed to eradicate order 
in typical diluted magnets (e.g. KMni_ 5 MgaF 3 ) [28,29]; 
in typical magnetic systems, the decrease in the ordering 
temperature and its saturation arc governed by the perco- 
lation threshold (where the ordering temperature vanishes 
at the critical dopant concentration of 5 C = 0.69 for the 
simple cubic lattice) . The faster degradation of orbital or- 
der with doping vis a vis percolation physics can indeed 
be attributed [11] to the directional character of the or- 
bital exchange interactions. A loss of orbital order was 
also seen in the bilayered ruthenate Ca3Ru207 [12] with a 
suggested critical pressure of P ~ 55 kbar; this degrada- 
tion of orbital order correlates with a metal to insulator 
transition. Our work demonstrates the possibility that 
quantum criticality is associated with the vanishing orbital 
ordering temperature of systems such as KCui_5ZnaF 3 
and may be accessed, and its signatures probed for, by 
varying parameters such as the external pressure and dop- 
ing. Orbital-only interactions (cither borne directly by the 
orbital-only Jahn- Teller effect or by the orbital component 
of the orbital-dependent spin-exchange) are pertinent to 
numerous other systems. These include, for example, the 
layered cuprates (K2CUF4), and manganites (La2Mn04 or 
La 2 Mn 2 07). Such systems were discussed recently in [23]. 

Note added in proof: Several months after the initial 
appearance of our work, Ref. [23] discussed non-uniform 
structures that may appear in doped orbital systems de- 
scribed by orbital models. 
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